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Fundamental superstrings (F-strings) and D-strings may be produced at high tem- 
perature in the early Universe. Assuming that, we investigate if any of the instabil- 
ities present in systems of strings and branes can give rise to a phenomenologically 
interesting production of gravitons. We focus on D-strings and find that D-string 
recombination is a far too weak process for both astrophysical and cosmological 
sources. On the other hand if D-strings annihilate they mostly produce massive 
closed string remnants and a characteristic spectrum of gravitational modes is pro- 
duced by the remnant decay, which may be phenomenologically interesting in the 
case these gravitational modes are massive and stable. 

PACS numbers: 11.25.-w,11.25.Wx,98.80.Cq 



I. INTRODUCTION 

After the dismissal of cosmic string as candidate seeds for the formation of structure, 
they received a fair share of attention in the last few years from both the theoretical and 
observational perspective. 

Scenarios of inflation involving branes have become popular recently , because of the 
appealing feature of providing a working model of inflation based on a fundamental theory 
on one side, and because they can represent a phenomenological handle for string theory on 
the other. In this kind of scenarios the production of D-strings as cosmic defects, by Kibble 



mechanism, is a generic prediction |3j], even if their (meta-)stability is not at all ensured 

In general, cosmic defects like monopoles, strings and walls can be problematic as their 
red-shift is slower than the one of radiation: if present in the early Universe and if cosmo- 
logically stable they tend to overdose the Universe, unless they are diluted by some release 
of entropy into the Universe after their creation. While this is certainly true for stable 
monopoles, whose annihilation rate is negligible for standard choice of parameters Q], in 
the standard scenario cosmic strings loose energy by chopping off small loops and finally 
reach scaling solutions on which they make up a constant fraction of the energy density of 
the Universe during both the era of radiation and matter domination This allows the 
density of long strings to be always roughly few per Hubble volume, while the small loops 
disappears by emitting gravitational radiation. 

The situation is less clear for domain walls, as the existence of scaling solutions and the time 
required for the solution to fall into these kind of attractors in presently under investigation 

0. 

In this work we concentrate on cosmic strings. In particular we want to identify cosmic 
strings with the D-strings of string theory. 

The production of lower dimensional defect can result for instance from Dp Dp-brane 
annihilation, which will produce D(p — 2n)-branes as remnants. In a type II theory Dp- 
branes are BPS only if p has a definite parity (even in IIA, odd in IIB), Dp-branes with the 
"wrong" dimensionality will decay immediately (on time scale given by the string scale). 

Actually in type I string theory non BPS DO-branes can be stable, because the unstable 
mode is projected out by the orientation projection taking type IIB into type I. Still it 
has to be remember that the presence of a type I excitation spectrum is related to the 
presence of orientifolds. Even if the starting setup is endowed with a space filling orientifold, 
since some of the dimensions may get small during cosmological evolution, T-dualities over 
those directions getting smaller than the string length confine the orientifold hyperplane to 
subspace of strictly positive codimension, so the bulk physics, away from the orientifolds, will 
be again unoriented, i.e. of type IIA or IIB. Thus even taking a full type I configuration as 
a starting point, the general assumption that only the BPS Dp-branes are stable is justified 
on general grounds. 

Strictly speaking BPS D-strings are present in type IIB and type I theories and not in type 
IIA, but we can consider Dp-branes with arbitrary p wrapped around (p — 1) dimensional 
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cycles to be D-strings from the 4 dimensional point of view, allowing a variety of values for the 
tension of the effective D-string. For instance, the tension of a Dp-brane 
(being a' = l 2 s the square of the string fundamental length) and assuming that a D(2 + n)- 
brane wraps around a n dimensional flat torus with radii R\ , . . . , R n the effective tension /i 
of the resulting 4 dimensional D-string is 

1 ^rRi 



2ii l s g s fjr l s 



n 



i=i 

Moreover the internal dimensions do not need to be flat, but can be warped (a la Randall- 
Sundrum jjg:], for instance), so that their tension can be redshifted (or blue-shifted) by huge 
factors and in principle it can take any value. 

One important and not completely settled issue is the stability one. The (effective) D- 
strings are generically created but for them to play a role in cosmology they must be stable, 
or at least metastable over cosmological time scale. 

n 

An instability already pointed out in |9] is the following. D-strings couple electrically 
to a 2-form, whose dual in 4 dimensions is a scalar, axion-like field. Once supersymmetry 
is broken the axion acquires a periodic potential with several degenerate minima, so the 
axion is expected to condense to different minima in different regions of space, separated by 
domain walls. Finite domain walls are bound by D-strings, and the domain wall tension, 
unless it is exceedingly small, will make D-strings shrink and collapse. 

Another source of instability is breakage, which is present whenever systems of Dp and 
Dg-branes are present at the same time. Let us define v as the number of dimensions in 
which the Dp branes are extended and the Dg are not or vice-versa, i.e. the total number 
of Neumann-Dirichlet (ND) plus DN dimensions. If v = 2 the axion symmetry is Higgsed, 
there are no axion domain walls, so the branes can break. Strings can in general break or 
annihilate if they meet a brane with v = 2, if v = 4 or v = 8 the system is stable and 
super symmetric, but we fall back into the domain wall problem when supersymmetry is 
broken (if v = 6 branes repel each others). 

Since the branes come with orientation, meeting a brane with opposite orientation they 
can also annihilate, whereas if they have exactly the same orientation and u = 0,4,8 the 
full system is still BPS. Of course intermediate situation are possible, when branes meet 
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e. A description of the microscopic mechanism leading to recombination is given in 
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The stability issue is thus far from settled but it can be solved in specific models. Here 
we will focus on the effects triggered by two of the aforementioned sources of instability, the 
meeting of a D-string at angles and the annihilation of D-strings. This kind of processes 
can be addressed in the context of string theory and the dynamic of the instability can be 
studied through the time evolution of a specific tachyonic mode of the system 
Since everything is coupled to gravity such instabilities can be expected to give rise to some 
gravitational radiation, which we study quantitatively in the present work. Modeling the 
string sources with an effective energy momentum tensor the amount of emitted gravitational 
radiation can be computed, but a negligible result is found. 

Actually, rather than the emission by a time dependent dynamics we can consider the 
closed string remnant left by a brane-anti brane annihilation Q], whose decay can be an 
interesting source of gravitational radiation, as we will show. 

It is also worth noticing that no cosmic string has ever been observed directly and that 
the cosmic microwave background constrains cosmic string tension /i to be GV/i < 1CT 7 , 
where Gat is the Newton gravitational constant. 

The outline of the paper is the following. In sect. |H] the salient features of the analysis of 
[lol | of the D string reconnection are reproposed and the GW emission in such a process is 
computed, then generalised to a cosmological setup like in Q|. In sect. IHII the annihilation 
process between a D string and a D string is analysed, finding again a negligible gravitational 
radiation emission. In sect. IIII Bl we study the gravitational radiation emitted by the decay 
of the massive closed strings which constitute the remnant of the D D-brane annihilation 
and conclusions are finally drawn in sec. IIV1 

II. GRAVITATIONAL RADIATION FROM STRING-STRING 

RECONNECTION 

The reconnection process of D-strings meeting at an angle and zero relative velocity has 



been studied in 10]. Adding a relative velocity between the strings boils down to introducing 
a non one (and non zero) probability of interaction jlll . Il^j ] . 

Stacks of parallel D-branes admit a low energy description in terms of a U(N) super Yang 
Mills theory, where iV is the number of D-branes into play. The bosonic excitations of the 
F-strings stretched between the D-strings are gauge vectors (scalars) from the point of view 





Figure 1: Two D-branes at an angle 7 before and after recombination. 

of the D-string Poincare group, if polarised along (orthogonally to) the brane. The system is 
described by a 1+1 dimensional Yang Mills theory coupled to as many scalars as transverse 
direction. The scalar fields parametrise the displacements of the brane in the orthogonal 
directions. 

A. Almost parallel D-strings 

1. Reconnection dynamics 



Following the analysis of [lfj , let us consider a system made of 2 almost parallel D-strings 
at an angle 7. The low energy dynamics of the system is described by the action 



S = -| (2W) 2 Tr J dxdt Q V" + -D&D^ 



(1) 



where \i — 0, 1 and i E {2 . . . d} runs over the dimension orthogonal to t and x. 

In order to describe two D-strings at an angle 7 the previous action has to be expanded 
around the background $>f' A$ 



$f > = 6 h2 qx a 3 , Af = 



(2) 



with q = tan 7/(^0!') and a % the Pauli matrices which form the algebra of 577(2), plus the 
fluctuations $ (/) , A^ 



(/) 



fc,2 {Z(t, x) a 3 + <f)(t, x) a 1 ) , A<p = 6^1 a(t, x) a 5 



(3) 



To lowest non trivial order in interactions, the resulting Lagrangian density is, |31| 

c ~ id 2 + 1 (> + - 1 (y 2 + r 2 ) - y - ?e ; - ^ +0+^(9+ o 

-ia 2 ( g V + 2 5 a£) , 
and it takes to the equations of motion 



(4) 



a + qxcj)' — q<f) + aq 2 x 2 = , 
4> — <fi" — a'qx — 2aq = , (5) 
£ - C" + «'0 + 2a 0' + a 2 <?z = , 

where it is understood that a dot stands for a time derivative and a prime for a derivative 
with respect to x. 

We note that consistency of the equations requires that a, <p are of the same order, while £ 
is 0(a 2 ), so that the first two equations of © can be solved for a and <p and then the result 
plugged into the equation for £. This justifies a posteriori the choice made in (0J to drop 
0(a 3 ) terms as well as 0(a£ 2 ). 

For the coupled a, system it is convenient to take the ansatze 

a(t, x) = J2 n>0 A n e imnt a n (x) , 
0(t,x)=En>o^e im "*a n (x), 

where the equation of motions force the condition m 2 = (2n — l)q. 

Both a and have an unstable lowest lying mode, with a Gaussian profile, so that neglecting 
the positive mass squared mode the solution can be written as jlO[ 

a(t, x) = A e^ qt e- qx2/2 = 0(t, x) . (7) 

The expression for £ is slightly more complicated, but again its non trivial part is concen- 
trated in an area of size q' 1 ! 2 around x = 

£(t,x) = A 2 Q e 2 ^*/^ [e 2 ^ x erf(^x + 1) + e - 2 ^ x erf(^x - 1) - 2 sinh(2^x)] . (8) 
2 y q 

In the previous (JSJ) the integration constant of the solution to the homogenous part of the £ 
equation is fixed by requiring that the lim^-too £ = independently of t. 

To give a geometrical interpretation we remind that the $j's keep trace of the displace- 
ments of the D-strings in the orthogonal directions, their eigenvalues are proportional to 
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the locations of the D-branes in the z-th direction. In our case the distance d y along the y 
direction between the 2 D-strings is d y = 2 x 2i\ct^ y where <& y is given by 



$= l/gx + e 0o(x) \ 1 / ^q 2 x 2 + 2qx£ + 4> 2 ^ 

2 \ 0q(x) — gx — w \ ■>/ q 2 x 2 + 2gx£ + </>q 

where in the previous a gauge transformation has been performed to diagonalise the 
field and only the first order in £ has been kept. The dynamics of the Yang Mills system tells 
the D-strings to recombine smoothly around the intersection point. This analysis neglects 
in the Lagrangian terms of the order a 4 and higher, thus it is valid for 

A e^~ qt e- qx2/2 < y/q. (10) 

By a shift in time we can set Aq = y/q (as we will do in the rest of the paper), so the solution 
(J7J) can be expected not to be valid for t > 0. 

Here we neglected the effect of a relative velocity v between the D-strings, which affects the 
probability of the reconnection taking place. It has been studied in with the result that 
the probability P ~ 1 if v <C 1, which takes to the condition 



v ^e^«Vg. (11) 
Again this restriction suggests that this picture is a good description only for negative times. 



2. Energy momentum tensor 

Having understood the dynamics of the recombination via the effective action of the 
lightest string modes, it is then possible to compute the energy momentum tensor of the 
system to estimate the amount of gravitational radiation emitted during the process. We 
take the definition 

T^ = jL— . (12) 

Explicitly on the solutions (JZJ) the t, x components of the energy momentum tensor are 



T/jj, = jj,2ir 2 a 2 



-<7 2 /2 - q? -qi 

-qi a 2 q - g£' - q 2 /2 J (13) 
x [S (y — qxna) + 5 (y + qxTta')} S(z) 
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where it is understood that the strings are immersed in 4-dimensional space time. 
In eq. (jl3j) we recognise the constant contribution from the background configuration of the 
tilted D-strings, which is T-dual to a constant magnetic field on a D2-brane. The time 
dependent part is the source of the gravitational radiation. This energy momentum tensor 
is conserved, but it would not have been conserved without the inclusion of £ in the ansatze 
(J3J), like in [10| . where the solution considered have £ = 0. This should not come as a surprise 
as forcing the solution © with £ = would imply an explicit breaking of Lorentz invariance, 
which would necessarily let no right to expect a conserved energy-momentum tensor. The 
inclusion of £ turns the Lorentz symmetry breaking of the background into a spontaneous 
one, thus allowing the conservation of the energy momentum tensor. 
By linearising general relativity over a flat background 

g,*, = Viw + , (14) 

we have 

Dh^ = IQiiGnT^ , (15) 

where h^ v = h^ u — Xjlri^h^ has been defined and the condition dji^ = has been imposed, 
consistently with the conservation of the energy momentum tensor. Following 3| the energy 
momentum tensor can be Fourier transformed both in space and time 

f^(cu,k) = J T^(t,x)e i ^ t - l ^dtd 3 x, (16) 

to obtain the formal solution to (|15|) as 

h, u (t,x) = 4G N [ d*x'-J— f ^Mf^ufie-*^*-*^ 



2G 



N 



2tt(2*)' ^ ' ' (17) 
f„Ju,ujh)e iuj{t - x) du, 

TTX " 

where h is the unit vector pointing into the direction of x and the local wave zone approxi- 
mation \x — x*\ ~ x — x' ■ h has been used, with x = 
Given (fT3j) and the explicit solution (JJJ) 

a(t, x) = ^de^ qt e- qx2/2 for t < (18) 



one can compute tin, say, to be 

h n (t,x) = -GV/ig^W/ 2 f duo e -M*-*) 

x J 2^/q + iu l j ( 19 ) 

-G Ni i q^a'^e^^m - x, 9, 0) . 
x 



2 2 
LP r u; 

e 4< j + + e 4 « 
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valid for t < x, where f(t — x, 9, 0) takes account of the angular dependence of the radiation 
and its explicit form is 



f{u,6,(f)) = e 



+ e s - 



1 erf + 



(20) 



where 5± = | sin #(cos 0±g7ra' sin <ft) | and 6, <ft are the usual polar angles. Actually the details 
of the angular distribution are of little importance in view of the intensity of the radiation, 
as we will now show. 

To compare with experimental sensitivity it is convenient to express the previous result 
in terms of the (double sided) spectral density Sh(f) which for a short burst of (Fourier 
transformed) amplitude h(f) and duration At reduces to 

S h (f)~h 2 (f)/At. (21) 

From the computation (fT9*j) we can estimate a spectral density function, neglecting the 
angular dependence in the exponential, 



3 (GV/i)V/V 4 e - £ 



x 2 l+cu 2 /Uq) 
1n -i20tT-i/ + N 5/2 f G N fi Y ( a' \ V2 f x \ 



(22) 

2 / / \ 3/2 / \ -2 v > 



where redshift factors has been neglected. The typical signal frequency is ~ y/q/(n(l + z)). 
For presently working experiments Sh ^ 10 _42 Hz _1 [17], so the result (}22|) is ridiculously 
small. This is due to the inherently microscopic origin of the dynamics, whereas to generate 
a consistent amount of radiation a coherent motion of a macroscopic mass is needed. Let 
us now investigate if a cosmological population of such recombining D-strings may help in 
enhancing the signal. 



3. Cosmological rate 

One can expect that several processes like the one described in the previous section take 
place in the Universe and that the signal resulting from the sum of all of them is highly 
enhanced with respect to the one found in (}2*2*j) . 

We assume the Universe evolution is the standard one after cosmic string production, 
i.e. radiation domination followed by matter domination. Very sketchily, simulations of 
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the evolution of a network of strings show that no matter what is the initial density of 
strings, within few Hubble times they reach a scaling solution characterised by few long open 
strings per Hubble volume whose length stretches with the horizon, and a large number of 
small closed strings A network of long open strings has a negative effective pressure 
p p = —p p/3, thus it redshifts more slowly than both radiation and non-relativistic matter, 
but it does not overdose the Universe as open strings loose energy by chopping off small 
loops which eventually decay by emission of gravitational radiation. At an effective level one 
can write down a continuity equation which takes into account both the Universe expansion 
and the energy lost into loops 

Pop = -3H(p op + p op ) - = -1H pop - , (23) 

where < A < 1 is a number taking account of the possibility of a non-unity value of the 
probability of reconnection. For a generic Hubble parameter H = (3/t and p oc p/L 2 eq. ()23|) 
admits the scaling solution 

L= w^f- (24) 

allowing the string network to track the dominant energy source of the Universe and to 
make up a constant fraction of it. 

String loops have typical size l c \ and density n d roughly given by 

l d ~ at , n d ~ a'h' 3 . (25) 



The value of a is not known, as an indication we can take the value given in |19j a ~ TGnP 
where T ~ 50, so that it can well be that a < 1 (we remind that smoothness of cosmic 
microwave background requires GnP < 10~ 7 ). A relation similar to (|25|) between l op ,n p 
and t can be expected to hold for long, open strings with a ~ 1. Given (J25|) the number of 
reconnection events between open strings per unit of spacetime volume can be estimated to 
be 

Vop-o P (t) = n 2 op va p-op - Pt~ 4 , (26) 

where the cross section a op - p has been estimated as cr op _ op ~ P/L, with P the reconnection 
probability, see (jlljl . and relativistic velocities has been considered. Had we considered 
instead the interaction between open and closed strings we would have obtained 

v p-d{t) = n pn c ivo p-ci - Pt~ 4 a~ 1/3 , (27) 
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which is slightly higher than (|26|). where it has been used that the average distance traveled 

— 1/3 

by a string before meeting another one is ~ n 



'cl 



We want now to relate the time variable t to the redshift factor (1 + z) as in |14|. Let us 
consider the usual FRW metric 

rfs 2 = _ dt 2 + fl2(t) ^2 + r 2 d ^ (2g) 

the comoving distance r between an object emitting light at a time t(z), where (1 + z) 
denotes the redshift factor, and the observer at present time to is given by 

rW =£$>- (29) 

Parameterising the scale factor as 

1 /2 

<^e<j (t/teq) tj < t < t e q , (30) 

a (t/t ) 2/3 t eq <t<t 0} 

where £j mark the epoch of the beginning of radiation domination and t eq denotes the time 
of matter-radiation equality, the redshift-comoving distance relation is obtained 



3t 



;i + ze q y 1/2 



r K z ) = — x \ n . ^V2 



1/2 -, , 2 g — £ eg 



1 + ^-1 + - 



3 1 + 2 



2!j > Z > Z e q 



Yj2 Ze 1 



(31) 



(1 + z 

which, neglecting numerical factors, can be well interpolated by 

r (^)~^^_. (32) 
W a l + z v 1 

The previous formulae allow us to write the volume element as 

dV(z) = W(t )r 2 dr ^ 10 2 ^ (1+ ^ )13/2(i + ^ )1/2 ^ . (33) 

Let us introduce N(z), the reconnection rate for events occurring at a given redshift z. Its 
differential is given by 

dv m 2 n 4- ziz ) 3/2 

dN=(l + zT x v{z)z^-d\Kz ~ -z 3 1 + ; 1 rflnz, (34) 

^ t (1 + 2) 3/2 

where the time-redshift relation 

t(z)=t (l + zy 3 / 2 (l + z/z eq r 1 / 2 (35) 
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has been used. The observer at r = sees a burst with duration At = At(l + z), being At 
the proper time burst duration. By comparing (At ) _1 with dN/dlnz we can verify if the 
signals actually overlap to make a continuous background, or if they are resolved in time. 
Since N(z) is a growing function of z, it turns out that for z* > z the signal is made by 
a collection of bursts whereas signals originated at higher redshift than combine into a 
stochastic background. Comparing (J3*lj) with At one obtains 

^~10 10 (tan 7) 1/8 a 1 ' 4 

The average over 7 introduces a numerical factor ~ 0.9. In the case of continuous back- 
ground, i.e. for events generated at redshift z > z*, taking into account the cosmological 
redshift, the spectral density function can be estimated as 

S h (f) = £ n \h(m + z'))\"^-d\nz' , (37) 

where z in stands for the red-shift of the set in of the scaling solution and h(f(l + z)) stands 
for a sum over metric polarisation and average over directions h of h fiu (f(l + z),n). Since 
dN(z)/d\nz oc z 3 at large z, the integral in eq.()37jl can be approximated by the value of the 
integrand at the upper limit of integration region 

in 2 z 3 

S h (f) = \Hf(l + z in ))\ 2 ^^, (38) 

*o z% 2 

where h(f) can be read from (|2*2*|) with the aid of (}2*T]) after averaging over orientations. 
With respect to the signal from an astrophysical source given by (|22|) there is a factor 
10 2 zf n / (to z z Jq) instead of At — g~ 1 ^ 2 (l + z), leading to an enhancement for high enough 
redshift, but again the final result is negligible 

AC/) = 10-^-^4 (^tf (^) 2 (^) 3 , (39) 

where an average over D-string meeting angle 7 has been performed and the distance x 
from the source has been expressed as a function of the red-shift according to (|32|1 . Moving 
to higher redshift brings competing effects: a higher distance from the source, see eq. ()32|) 
and a decreased integration volume, see eq.(|HHJ), both of which tend to suppress the signal, 
and increased event rate, see eq.([26|). which tends to increase the signal and it is actually 
the dominant effect. We remind the relation between redshift z and temperature T during 



a 



11/2 



TeV 



-1 



(36) 
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standard cosmological evolution 



Actually before the scaling solution sets in, the density can be higher than (J25)l and 
consequently the rate of string encounter higher than (]26[). One could consider the very 
extreme case of an initial density of cosmic D-strings 

n ex ~ /^V" 3 (41) 

at the epoch of cosmic string creation (inflation ending), characterised by a redshift z ex . 
Such a density would correspond to a Hagedorn energy density ph ~ a' -2 , which is generally 
a limiting energy density for open string systems. The corresponding rate of reconnection 
events per volume is 



(ex) _2 /-4 



op— op 



(42) 



while the standard scaling solution is recovered in a few Hubble times. The highly non 
standard form of the D-brane tachyon action, which suppresses gradients of the field, allows 
us to assume a far higher density of defects than the standard Kibble argument would allow 
as we will see in the next section, see [l^. With this assumption we can reconsider the 
integral (J37|) . noting that it will be dominated by the early epoch z' ~ z ex thus giving, 
analogously to (J3HJ), 

S h (f) ~ 10 3 V V ^ 2 ^ 1Q2 fe^ 2 ~ 10 ^ Hz -i e -*W4. f Vf -%>43) 

which is still way too small for detection. It has also to be noted that the Kibble density is 
smaller than the extremal density n ex given in (j4*Tj) only if 



< G N 1/3 l;^ . (44) 



With the process considered so far one cannot obtain a large enough result, as an in- 
herently microscopic mechanism is involved, which affects a portion of the string of size 
the order of the string length in a huge Universe. For the gravitational emission to give a 
sizeable effect we need it to be produced in a coherent motion of all its world volume. In 
sec JIIII we move to the study of systems of D-string D-string which annihilate, but before 
doing that we complete the analysis of recombination by describing two almost antiparallel 
D-strings, which is a system close to the D D-string case. 
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B. Almost antiparallel case 



1. System dynamics 



If the system of two D-strings is close to a brane-anti brane pair it can still be analysed 
through a Yang Mills low energy effective action, but a scalar tachyon has to be added to 
the spectrum, as it comes from the lowest lying excitations of the string stretched between 
the brane and the anti-brane Qj.The action is then 



<S = -iiVlnn') 1 I iltdx |Tr ( + X -D^D^ 1 + ^D^TD^T - ^-T 2 



(45) 



with m 2 = —1/(20;') . Now the expansion has to be taken around the background (again 
q = tan7/(7ro')) 



^ b) =6 h2 qxa 3 , A^ = 0, T® = , 



(46) 



and the fluctuations 



= , T« = r{t, x)a 2 



(47) 



In terms of the ansatze ()46|) and (}4Tj) the Lagrangian, up to lowest order interaction, can be 
written as 



r 1 -2 1/2 1 2 1/22,0 t\ 2 , 1 2 , lf2 1 t /2 

£ oc -r r o ox + 2qxf r H r + -t c 

2 2 2 H 2 vy y v 4a' 2 s 2 s 



(48) 



The analysis follows straightforwardly as in the previous case to obtain the equations of 
motion 



f-r" + q 2 x 2 r-r/(2a') = 
£-£" + tV = 



(49) 



Again the unstable mode solution is 



e \/it 



r(t, x) = C Q e pt e- qx2 / 2 , 
cosh(2/?x) - i (e 2/3x erf (1 + /3x) + e" 2/3:E erf (-1 + /3x)) 



(50) 



where /3 ee v /l/(2a') - q. 
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Figure 2: Two D-branes at an angle tt — 7 before and after recombination. 



2. Gravitational radiation 



From the dynamics exposed in the previous section the energy momentum tensor can be 
estimated to be 



It is then straightforward in full analogy with sec III A2l to estimate the gravitational radia- 
tion, obtaining the same result. 



We have seen that the recombination process of D-strings described in sec. |H] is not 
a phenomenologically interesting source of gravitational radiation. This was due to the 
smallness of the size of the interaction region, which was of the order of the string length. In 
general interesting sources of gravitational radiation are related to coherent motion of large 
and extended portions of mass. 

We thus now turn to study a process in which two D-strings meet and annihilate, thus 
enabling a much larger reservoir of energy to be converted into gravitational radiation. Such 
systems have been studied in detail in a full string theory context for different configurations 
of the tachyon describing the instability jl^ . 




(51) 



x [S (y — qxna) + 8 (y + qxna')} . 



III. STRING- ANTI STRING ANNIHILATION 
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In particular it has been shown that the tachyon can acquire different spatial profiles, trig- 
gering the decay of the original system into different final states. For instance a trivial 
tachyon profile in the space dimensions will lead the brane anti-brane system towards the 
so called tachyon matter which will be discussed in sec JIII B( whereas a kink profile will let 
a Dp Dp-brane system decay into a D(p — l)-brane, a vortex configuration into a D(p — 2)- 
brane and so on. We show next that in the case of a trivial spatial configuration, important 
conclusions can be drawn for the production of gravitational radiation. 



A. Recombination and annihilation 



For ordinary phase transitions the presence of defects is related to the global configuration 
of some field triggering the transition. In an expanding Universe with Hubble rate H, the 
defect separation d must be smaller than the Hubble length H" 1 , as any super- horizon 
correlation would violate causality, thus leading to the lower bound of one defect per Hubble 
volume (Hubble density), and the actual number of defects per Hubble volume is as large as 
a few. 

However the density of defects can be much larger than the Hubble one as this picture 
does not apply to the tachyon condensation, whose dynamics is described by the Lagrangian 

s = -vf ^^(Tiyr+Wr, (52) 

where the V(T) is a symmetric, positive potential, which vanishes asymptotically for T — > 
±oo as e~ T |2jj fit is not very important here the exact form of such potential). 

As shown in 18], the non-standard form of the kinetic terms implies that once the tachyon 
has rolled to large values, gradients are exponentially suppressed, essentially eliminating the 
restoring force which would tend to eliminate gradients within a casual volume, thus allowing 
a higher defect density than the Kibble value. 

We assume that the density n ex as in (J4T|) can be achieved, corresponding to an energy 
density of the order of the Hagedorn scale, which is a limiting energy density for open 
strings Q. 

Actually in the evolution of an ordinary string network, the initial conditions are not 
crucial as the network reaches the scaling solution in a few Hubble time, long strings loose 
energy by chopping off small loops, which in standard simulations then looses energy by 
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emitting gravitational wave until disappearance. If we take the extremal case of Hagedorn 
energy density like in the end of sec. (jll A3|) the long string inter-separation is d ex ~ p^^a 1 
rather than the Kibble value dx ~ H" 1 ~ G~^~' 2 pT 1 , which is eventually reached after few 
Hubble times. 

The high starting density of strings is anyway not enough in itself to produce a large 
amount of gravitational wave if the usual picture is assumed, i.e. that small loop radiates as 
in ^| via the formation of cusps and kinks. Again this is due to the smallness of the physical 
size of the radiating objects, in this case the small loops, whose typical size is some orders of 
magnitude smaller than the Hubble scale, see eq. (f23|) . where H ~ y/G^p — lpi/l 2 , taking 
for p the extremal Hagedorn density. [32]. In this picture the emission is originated in the 
early Universe, when the size of cosmic string is small as they had no time to stretch with 
cosmological expansion, so a small amount of radiation is produced. In [14] for instance, a 
typical signal close to experimental sensitivity at a frequency /, with rate N, originated at 
redshift z is obtained from a kink in a cosmic loop of size l c i roughly given by 

l d * 10 20 m (1 + z)^ 2 (1 + z/z eq )^ 2 (^) , (53) 

and the rate A" of such an event at a given frequency / is roughly given by 

N = 2x 10 - V a r - ^) + "^"" (A) ^ (4-) ' m ■ (54) 

y (i + ^) 7 /6 ViO" 6 / v kHz / 

B. Annihilation with spatially trivial tachyon profile 

So far we have dealt with an effective field theory description of D-string D-string an- 
nihilation, the excitations of the branes have been described as classical fields subject to 
ordinary second order differential equations. 

On the other hand the annihilation process admits an exact description in terms of the 
fundamental string theory, i.e. in terms of the two-dimensional world-sheet theory, by adding 
to the standard world sheet action 



S ws = -^J drda (^-d a X»d a X»g^ + Vp a d a r9^ , (55) 
a boundary term 

^ocA / dr^°e x °^ ® a x . (56) 
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The world-sheet (disk) partition function Z ws , with the world sheet couplings interpreted as 
space-time fields, corresponds to the space time action S 23], given by 

S = Z ws oc J cFxy/^g J [dX") [dip"] e - Sws ~ Sb . (57) 

The boundary term describes the condensation of the tachyon to a profile of the type 

T = Ae*/^W ® <r x . (58) 

Giving an expectation value to the tachyon shifts the mass squared of the open string modes 
by an amount Am 2 oc (T 2 ), eventually dropping them out of the physical spectrum as the 
tachyon increasing expectation value makes them more and more massive. 

A D-brane in the free theory can be described by a specific boundary state \B), the 
\B) — > | B) vacuum amplitude is written as the so called cylinder amplitude 

A BB = (B\P\B), (59) 

where P is the closed string propagator. When the internal particle running in the propaga- 
tor goes on shell the amplitude Abb picks an imaginary part. Once written the expansion 
of \B) into a set of complete states as 

|3) = J>("/) |/> , (60) 
/ 

by the optical theorem we have 

lmA BB = J2^\U(u Jf )\ 2 , (61) 

where u)f is the energy of the state /. In the unperturbed theory, where the tachyon has 
vanishing expectation value, the brane boundary state can be written as 

\B) = 0\Q,k), (62) 

where O is an operator build out of fermionic and bosonic matter oscillators and ghost 
oscillators and \Q,k) is the Fock vacuum with momentum k. The specific form of the 

nn 

vacuum will not be needed here, as we shall simply quote known results, see [20, 21]. Taking 
the initial condition to be that of D-strings, the addition of the boundary term (|56|) turns 
\B) into \B) T , whose energy momentum tensor is obtained from (|H7j) via ()12j) 

7V = 2/idiag(-l,£(t),0,...,0) , (63) 
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with 

9 ® ~ l + 2 7t 2 \ 2 e^ t / l ° ' ^ 

The boundary state \B) T at t — > — oo corresponds to \B), and at t — > oo tends to a state 
without the original D-strings but endowed with the same amount of energy and vanishing 
pressure. 

Applying the optical theorem as in (JfiTjl one can find, see Q j that the expansion 
coefficient Ux{oJf) relative to \B) T turns out to be just the Fourier transform of g(t), which 



is 



PtM\ = . J 1 ^ ■ (65) 



The imaginary part of y(5 |-B) T is related to the average numerical density n and energy 
density p of particles emitted by 



EttH 17 ^/)! 2 ' ( 66 ) 



n 

/ ; 

/ 

where the sums run over a basis of closed string states, implicitly including both the sum over 
excitation level iV of the closed string and the momenta of the spatial directions transverse 
to the brane. For large iV the degeneracy of states Q(N) of a closed string in D dimensions 
is 

0(N) = N -(D+D/2 e 4*y/Z=SVN ^ (6g) 

and here we assume D = 10 from superstring theory. Given the degeneracy of the states 
and the mass formula for closed strings 

/ 2 ■ (- 1 



aW = 4 \N - - j , (69) 

the exponential in eq. (|63j) and (p3%|) combine to leave a power law spectrum for the radiated 
energy according to 

p oc / dE E- p/2 , (70) 
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for Dp-branes, with a total emitted energy which is ultraviolet divergent for p < 2. Of course 
the emitted energy cannot be larger then the mass of the initial brane, so back-reaction is 
expected to set in and cut off the emission, with the result that the total radiated energy is 
finite and of the order of the initial mass of the brane. In particular cutting off the integral 
in eq. (f7T)|) at E cut ~ pl ov would give all the energy into closed string modes, emitting most 
of their energy into very massive strings. 

Actually eq. (fTUj) seems to suggest that for p > 2 the amount of energy of the initial Dp- 
brane going into closed string modes is rather small. However the decay described by the 
boundary term (jHEJ) is spatially homogeneous and it is not expected to realistically describe 
the decay process. Perturbations with a non trivial spatial profile can still be tachyonic and 
thus contribute the decay. We should expect the decay to take place rather inhomogeneously 
and the original Dp-brane can be thought of decomposing into a collection of small patches. 
In each of this small patches, disconnected one from the other, the behaviour of a Dp-brane 
is similar to that of a collection of DO-branes 1221. 



C. Gravitational mode production 

The problem of gravitational emission has now been traded from the study of annihilating 
D-branes to the one of decaying of a highly massive (hence highly excited) closed fundamental 
string. Two-body decays of F-strings have been studied in literature for both the bosonic 



string |26[ and the superstring 27, 28]. The decay of a closed string into n closed strings is 
a perturbative process whose amplitude r oc g™, thus in the perturbative regime g c <C 1 the 
dominant process is two-body decay. 

Let us consider the amplitude decay A of a specific initial closed string state 4>n of mass 
M ~ 2\/~N /l s into two given final states 0at i; 4>n 2 of masses respectively mi >2 ~ 2^J Ni^/l s - 
The inclusive decay width of a generic fundamental string excited at mass level N, into two 
generic final states at mass level Ni,N 2 is obtained by integrating the amplitude squared 
over the available final states, and averaging over initial states at mass level N. We assume 
to have d — 3 + 1 large ordinary dimensions and d c small ones, and a highly energetic (larger 
than the string scale) closed fundamental string which can both wrap and have momentum 
around the small dimensions. Phase space integration includes a sum over winding and 
momenta. 
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The differential decay width dT d ^ mi 2 is then given by averaging over initial states and 
summing over final ones 27 1 

« = mm) £ w ,<w *-* K w (jV - Ni - eJ Wr L d * h -" • (71) 

The two-body phase space in d = 3 + 1 dimensions is 



d ^+P = -J^—dn (72) 

h,h 16tt 2 M 



where is the modulus of the momentum of the decay states and Q the solid angle. Con- 
sidering for simplicity flat internal compact dimensions, the inclusion of momentum and 
winding modes in the mass formula (|6T)|) adds a sum over momenta and windings in the 
compact directions. The closed string mass formula which generalises (joTJjl in the presence 
of discrete momenta and windings is 

aW = 4 (iV-0 + ^( a '| + ^ , (73 ) 
which for large energies can be written as 



This allows to approximate the summations over rii and Wi as Gaussian integrals in the case 
respectively m ^> hR^ 2 and m 3> R 2 /l^ s - Using the asymptotic formula for the density of 
states (|68|) the amplitude (|7Tj) can be integrated to 



2 (2n/a) d ° hE /mjm^-^-W) 2 ( 2E mi m 2 \ 



d-3 

2 



(75) 



TrM 2 ~ dc l 3+D ~ d ~ 2dc 

where a = Air/y/S, d c the number of compact dimensions along which the closed strings have 
momenta and windings, the kinetic energy of the decay objects 

E = M -mi-m 2 ~ — — — k 2 (76) 

2mim,2 

has been introduced and the approximate relation 



N - Ni - N 2 ~ a'M\/k 2 + m\ - m\ ~ am x m 2 (77) 

has been used. 



22 



Let us call T m the total rate of production of a string of mass m, which is obtained by 
summing T over all values N 2 . The sum over N 2 can then be trade for a continuous integral 
over E with the substitution Y1n 2 ~^ ~a'm 2 dE, obtaining 

rS mi , <x ^ "-y W-" 72 / B*^-*'^ , (78) 

where tur = m(M — m)/M. Since in an interval AN corresponds to a'mAm/2, the density 
of states per mass level is p(m) = a'm/2, thus we can define the differential decay width 
per mass interval j^] 

= rE roi , 2 p(m) oc glM (l s Mf~ d {l s m R y^ . (79) 

The differential decay rate ()79|) clearly shows that light strings are more easily produced 
than heavy ones. Thus we can safely assume that a consistent fraction of the initial rest 
mass is converted into massless closed strings, i.e. radiation. 

We are then lead to consider the amplitude inclusive over m 2 and with mi = which is 

H 

dT cl _ rad ~ gl ^^J sMiUY^dk , (80) 

similar to a black body spectrum, leading to a total decay rate 

r d ^ rad ~g 2 c M. (81) 

Thus in a time much shorter than the string scale a D-string and a D-string decay efficiently 
into radiation. 

If it is so one of the two decay objects can still be massless from the ten-dimensional point 
of view but it will be massive from the effective 4-dimensional theory. In this case the decay 
width ()80|) should have been modified by the inclusion of an additional integral 

dV d ^ KK = gl-^—— m ^l s M{l s k) d - 1 dk n Rtdki ~ g 2 c e- ahVW +^m d ^ 1 dm . (82) 

^ i 

Then integration over the momenta in the compact dimensions have the effect of multi- 
plying the differential decay width by a factor Y[Ri/h, for Ri 3> l s , because of the opening 
up of new decay channels, and it does not modify the result (jHH|) for Ri ~ l s . 

If two D-strings are not exactly anti-parallel but meet at an angle 27, as in fig. El the 
tachyon condensation makes the D-string disappear not over the entire world-volume of the 
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L=21/tariY 



Figure 3: Recombination of two almost antiparallel D-strings. It is highlighted the region of space 
where tachyon condensation occurs. 

D-strings, but only as long as the D-strings are separated by a distance shorter than 21 s , 
which for strings meeting at angle 27 happens for a region of size L x = 2l s / tan 7, see fig. E3 
n he energy density in gravitational waves p gw is related to the spectral density function 

by ' 



GNPgw 



dffS h (f) . 



(83) 



The energy a single annihilation of D-branes make available for radiation can be estimated 
as \xL x which eventually expands spherically. From (j8*Uj) and (|83|) one can then compute the 
strain due to a gravitational wave emitted at a distance x from us (here we neglect red-shift 
for the moment, thus having fl s <C 1 at interesting frequencies) 

-al s f (Q A / \3/ 2/ 



_ ... 2 G N pL x a' l s f e~ alsf 
bhU) - - 
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al s f 



10 



-108 



Hz" r 



tan 7 VIO^ 7 J V T eV 



a 



x 
Mpc 



.(84) 



Of course Sh{f) cannot diverge for 7 — > 0, as the naive estimation for L x does, but for small 
7 one has to substitute 21 s j tan 7 with l op , the (open) cosmic string length. In this optimal 
case of perfect anti-alignment one has 

2G Nf il op a'l s f(l + z)e- al °^ ^ 

(85) 



st\f) 
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x- 



1 



10- 70 Hz- 1 e- a/s/(1+2) - 



,-al s f(l+z) 

a + z) 



1/2 



/ G N fi \ 



a 



r-2 



z*(l + z/z eq ) 1/2 VlO-V VTeV" 
where it has been used that l s f(l + -C 1, as a frequency I' 1 at an epoch characterised by 
a temperature T = Ij 1 , once red-shifted corresponds to 100GHz today, assuming standard 
adiabatic evolution in between. 
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Assuming that the radiation is produced by an annihilation event during cosmological 
evolution when the scaling solution has set in, the expected rate of such an event is 

N=^f(l + z) 3 / 2 (l + z/z eq ) 1/2 . (86) 

It is clear that here again we have to assume an extremal initial density of strings well 
above the Kibble density to hope to have a detectable signal. Let us then assume an initial 
density of annihilation ()42|) at an epoch characterised by a red-shift z ex . Annihilation will 
thus take place in a region of volume AV(z ex ), see eq. (1H3|) . and for a time duration 



At(z ex ) = 



dt 


Az ~ 


dt 






dz 


Z — Z ex 


dz 



to give 

Sf?(f) = ^- 2 G N d -^u ex AVAt = 
ft I dj 

2 G Nf il op a' I JQ- + z)e- al ^ 1+z ^ 1 t 4 z 4 

7T t 2 z 2 /{l + zf X 1 - e -«'./(i+*) X /xV 4 X (1 + z) 9 (l + z/z eq ) ~ ^ 7 > 

in 50 — Hz" 1 ( ^ ^ ( a ' \ e - alsS{1+z) 

(l+^)iV2(i +z /^)3/2 VlO- 7 G^V VTeV-V 

where we have considered the optimal case L x = l op . Since the extremal case of a Hagedorn 
energy density has been assumed, z ex has to be high enough to accommodate such a high 
density phase. Assuming standard adiabatic expansion after this extremal density phase the 
relation between the red-shift factor z and the string scale is eq.pU]) 

*= ioi5 (^r /2 - < 88 > 

which takes the estimate (|87|) very far from both experimental sensitivity and phenomeno- 
logical bound. In particular the nucleosynthesis bound requires that not too much radiation 
is present at the time of the synthesis of the primordial elements, in order not to spoil 
the beautiful agreement between theory and observation of primordial element abundances. 
Such bound constraints 

^(/)<10- 51 Hz- 1 (^ 3 . (89) 
Substituting (|55jl into (|57jl we can obtain 

s ' (/) = 10 Hz («f%) {m*} (10^5) 1 -e-^.. (90) 
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In our toy model the spectral strength can saturate the bound only around 104- 100GHz, i.e. 
at the end of the spectrum. This is hardly surprising as most of the energy of the F-string 
decay is concentrated for frequencies f' 1 ~ l s (l + z ex ). 

On the other hand if Kaluza-Klein modes are created and if they are absolutely stable, 
they can rapidly overdose the Universe. In this case of extremal production their initial 
energy density is a non negligible fraction of the total energy density, so their decay is 
necessary otherwise they would start dominating as soon as they become non-relativistic. 
Their final abundance is then set by the annihilation cross section j^^j], which may not be 
Planck scale, as gravity can be strong much below the Planck scale in large extra dimension 
scenarios. 

If instead we consider the standard scenario with the space-time rate of open string-open 
string encounters given by ([26)1. we can estimate the fraction of the critical energy density 
going into the production of massive gravity modes as follows. Let us suppose that each 
encounter give rise to the production of massive modes for a total energy Ekk = pL x , 
thus giving an energy density per Hubble volume pkk — pL x /t 3 . By normalizing it by the 
critical energy density p c = 3H 2 /(8ttGn) — (SttGa^ 2 ) -1 and integrating over the history of 
the Universe, using eq. (|3l)Jl . we have 

n KK (z nr ) ee P-BpA = 327 r 2 G N ^ r (1 + z'f 2 (1 + z'/z eq ) 1/2 dz 

Pc{Znr) to J z nr / n1 \ 

J \ 1/2 ( 91 ) 



a 



10 -5o4 ( Gnp \ 

z eq \TeV- 2 



■ 



where Zi n denotes as in eq. ()37|) the time of the onset of the scaling solution. Note that 
the integral in ()91|) has been halted at z nr , the red-shift at which the Kaluza-Klein modes 
become non relativistic. For z nr > z > z eq , p KK keeps earning a factor of z over p c thus 
resulting into a present fractional energy density Qrk = Qkk(z = 0) given by 

Z nr 1n -ll ( GnH\ I 01 



a "= 10-11 ^AsH ' (92) 

where for simplicity we have assumed the mass of the Kaluza-Klein modes to be a /_1//2 , see 
eq.(jH21), and substituted the relation (JHHJ) for z = z in = z nr . 

This result can be phenomenologically relevant for a still moderate value of the string scale. 
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IV. CONCLUSIONS 



Motivated by the recent revival in cosmic strings we have studied cosmic strings by mod- 
eling them as D-branes, the solitonic objects of string theory and computed the amount of 
gravitational waves emitted in recombination processes which can take place in the moder- 
ately early Universe. We have found that in general such processes are well below present 
sensitivities. Still the process leading to string annihilation can result into a production of 
gravity modes of fundamental strings which can be massive at an effective 4-D level and 
can thus have a cosmological effect. Such effect depends ultimately, in our simple model, 
on the string scale value, and it can be interesting even if the string scale is well below the 
Planck scale. From a more theoretical point of view, this mechanism represents a new way 
for producing out of equilibrium, very massive particles. 
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